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Abstract 

If  instead  of  the  full  motion  field,  we  consider  only  the  direction  of  the  motion  field  due  to 
a  rigid  motion,  what  can  we  say  about  the  three-dimensional  motion  information  contained 
in  it?  This  paper  provides  a  geometric  analysis  of  this  question  based  solely  on  the  constraint 
that  the  depth  of  the  surfaces  in  view  is  positive. 

It  is  shown  that,  considering  as  the  imaging  surface  the  whole  sphere,  independently  of 
the  scene  in  view,  two  different  rigid  motions  cannot  give  rise  to  the  same  directional  motion 
field.  If  we  restrict  the  image  to  half  of  a  sphere  (or  an  infinitely  large  image  plane)  two 
different  rigid  motions  with  instantaneous  translational  and  rotational  velocities  (ti,u;i)  and 
(t2,  W2)  cannot  give  rise  to  the  same  directional  motion  field  unless  the  plane  through  ti  and 
t2  is  perpendicular  to  the  plane  through  and  UJ2  (i-e.,  (ti  x  t2)  •  (wi  x  0:2)  =  0).  In 
addition,  in  order  to  give  practical  significance  to  these  uniqueness  results  for  the  case  of 
a  limited  field  of  view,  we  also  characterize  the  locations  on  the  image  where  the  motion 
vectors  due  to  the  different  motions  must  have  different  directions. 

If  (wi  X  0:2)  •  (fi  X  t2)  =  0  and  certain  additional  constraints  are  met,  then  the  two  rigid 
motions  could  produce  motion  fields  with  the  same  direction.  For  this  to  happen  the  depth 
of  each  corresponding  surface  has  to  be  within  a  certain  range,  defined  by  a  second  and  a 
third  order  surface.  Finally,  as  a  byproduct  of  the  analysis  it  is  shown  that  if  we  also  consider 
the  constraint  of  positive  depth  the  full  motion  field  on  a  half  sphere  uniquely  constrains  3D 
motion  independently  of  the  scene  in  view. 
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1  Introduction  and  Motivation 


The  basis  of  the  majority  of  visual  motion  studies  has  been  the  motion  field,  i.e.,  the  projec¬ 
tion  of  the  velocities  of  3D  scene  points  on  the  image.  Classical  results  on  the  uniqueness  of 
motion  fields  [6,  9,  10]  as  well  as  displacement  fields  [8,  12,  14]  have  formed  the  foundation 
of  most  research  on  rigid  motion  analysis  that  addressed  the  3D  motion  problem  by  first 
approximating  the  motion  field  through  the  optical  flow  and  then  interpreting  the  optical 
flow  to  obtain  3D  motion  and  structure  [2,  7,  13,  15]. 

The  difficulties  involved  in  the  estimation  of  optical  flow  have  recently  given  rise  to  a 
small  number  of  studies  considering  as  input  to  the  visual  motion  interpretation  process 
some  partial  optical  flow  information.  In  particular  the  projection  of  the  optical  flow  on 
the  gradient  direction,  the  so-called  normal  flow  [5,  11],  and  the  projections  of  the  flow  on 
different  directions  [1,  3]  have  been  utilized.  In  [3]  constraints  on  the  sign  of  the  projection 
of  the  flow  on  various  directions  were  presented.  These  constraints  on  the  sign  of  the  flow 
were  derived  using  only  the  rigid  motion  model,  with  the  only  constraint  on  the  scene  being 
that  the  depth  in  view  has  to  be  positive  at  every  point— the  so-called  “depth-positivity” 
constraint.  In  the  sequel  we  are  led  naturally  to  the  question  of  what  these  constraints, 
or  more  generally  any  constraint  on  the  sign  of  the  flow,  can  possibly  tell  us  about  three- 
dimensional  motion  and  the  structure  of  the  scene  in  view .  Thus  we  would  like  to  investigate 
the  amount  of  information  in  the  sign  of  the  projection  of  the  flow.  Since  knowing  the  sign 
of  the  projection  of  a  motion  vector  in  all  directions  is  equivalent  to  knowing  the  direction  of 
the  motion  vector,  our  question  amounts  to  studying  the  relationship  between  the  directions 
of  2D  motion  vectors  and  3D  rigid  motion. 


We  next  state  the  well-known  equation  for  rigid  motion  for  the  case  of  a  spherical  imaging 
surface.  We  describe  the  constraints  and  discuss  the  information  exploited  when  using  the 
full  flow  as  opposed  to  the  information  employed  when  using  only  the  direction  of  flow.  As 
will  be  shown,  whereas  full  flow  allows  for  derivation  of  the  direction  of  translation  and  the 
complete  rotation,  from  the  orientation  of  the  flow  only  the  direction  of  translation  and  the 
direction  of  rotation  can  be  obtained. 

The  2D  motion  field  on  the  imaging  surface  is  the  projection  of  the  3D  motion  field  of 
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the  scene  points  moving  relative  to  that  surface.  Suppose  the  observer  is  moving  rigidly 
with  instantaneous  translation  t  =  {U,V,W)  and  instantaneous  rotation  uj  —  (a,  ^,7)  (see 
Figure  1);  then  each  scene  point  R  =  (X,  Y",  Z)  measured  with  respect  to  a  coordinate  system 
OXY Z  fixed  to  the  camera  moves  relative  to  the  camera  with  velocity  R,  where 

R  =  — t  —  u;  X  R 


If  the  center  of  projection  is  at  the  origin  and  the  image  is  formed  on  a  sphere  with  radius 
1,  the  relationship  between  the  image  point  r  and  the  scene  point  R  under  perspective 
projection  is 


R 


Figure  1:  Image  formation  on  a  spherical  retina  under  perspective  projection. 


If  we  now  differentiate  r  with  respect  to  time  and  substitute  for  R,  we  obtain  the  following 
equation  for  f : 


r  =  ■ytr(r)  +  ■yrot(r)  =  — ((t  •  r)r  -  t)  -  u?  x  r  =  x  (t  x  r))  -  w  x  r 


Rl 


The  first  term  ntr(r)  corresponds  to  the  translational  component  which  depends  on  the  depth 
Z  =  |R|,  the  distance  of  R  to  the  center  of  projection.  The  direction  of  Utr(r)  is  along  great 
circles  (longitudes)  pointing  away  from  the  Focus  of  Expansion  (t)  and  towards  the  Focus  of 
Contraction  (— t).  The  second  term  nrot(r)  corresponds  to  the  rotational  component  which 
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is  independent  of  depth.  Its  direction  is  along  latitudes  around  the  axis  of  rotation  (coun¬ 
terclockwise  around  u?  and  clockwise  around  —oj).  See  Figure  2a,  b  and  c  for  translational, 
rotational,  and  general  motion  fields  on  the  sphere. 


(a)  (b)  (c) 


Figure  2:  Example  of  (a)  a  rotational,  (b)  a  translational,  (c)  a  general  motion  field  on  a 
sphere. 


As  can  be  seen,  without  additional  constraints  there  is  an  ambiguity  in  the  computation 
of  shape  and  translation.  It  is  not  possible  to  disentangle  the  effects  of  t  and  |R|,  and  thus 
we  can  only  derive  the  direction  of  translation.  If  all  we  have  is  the  direction  of  the  flow  we 
can  project  f  on  any  unit  vector  on  the  image  and  obtain  an  inequality  constraint: 


r  •  n,  = 


— ((t  •  r)r  —  t)  -  a;  X  r  1  •  n,-  >  0  or  <  0 


From  this  inequality  we  certainly  cannot  recover  the  magnitude  of  translation,  since  the 
optical  flow  already  does  not  allow  us  to  compute  it. 


In  addition  we  are  also  restricted  in  the  computation  of  the  rotational  parameters.  If 
we  multiply  a;  by  a  positive  constant,  leave  t  fixed,  but  multiply  by  the  same  positive 
constant,  the  sign  of  the  flow  is  not  affected.  Thus  from  the  direction  of  the  flow  we  can  at 
most  compute  the  axis  of  rotation  and,  as  discussed  before,  the  axis  of  translation.  Hereafter, 
for  the  sake  of  brevity,  we  will  refer  to  the  motion  field  also  as  the  flow  field  or  simply  flow, 
and  to  the  direction  of  the  motion  field  as  the  directional  flow  field  or  simply  directional 


2  Relationship  Between  the  Orientation  of  the  Flow  and  the  Depth-positivity 
Constraint 


If  we  have  the  flow  r,  we  know  the  value  of  the  projection  of  r  on  any  direction  and  we  set  all 
the  possible  information  by  choosing  two  directions  ni  and  112  (usually  orthogonal).  Thus 
we  have 

r  •  n,-  =  j^((^  ■  -  (^  X  r)  •  n,-  for  i  =  l,2  (1) 

We  can  solve  equation  (1)  for  the  depth, 

r  •  n,-  +  (oj  X  r)  •  n,-  1 


for 


z'  =  1,2 


((t  •  r)r  -  t)  •  n,' 

Knowing  the  value  in  both  directions  ni  and  n2  we  know  that  the  inverse  depth  has  to 

be  the  same,  and  also  has  to  be  positive;  thus 

r  •  ni  +  (u?  X  r)  •  ni  _  r  •  n2  -I-  (a>  x  r)  •  na 
((t-r)r-t)-ni  ((t  •  r)r  -  t)  •  n2  ^ 

If  on  the  other  hand  we  do  not  use  the  value  of  the  flow  but  only  its  direction  and  thus  the 
sign  of  the  projection  of  the  flow  on  n^-,  then  the  only  constraint  that  can  be  utilized  is  the 
inequality,  which  comes  from  the  fact  that  the  depth  is  positive.  Using  only  the  orientation 
of  the  flow  we  obtain  for  every  direction  Hj-: 


r  •  He  +  (u?  X  r)  •  n^- 
((t -rjr- t) -iii  ^ 

This  inequality  provides  inequality  constraints  on  the  rotational  and  translational  compo¬ 
nents,  which  are  independent  of  the  scene:  If  we  consider  the  sign  of  the  translational 
component  ((t  •  r)r  —  t)  •  n,-  and  the  sign  of  the  rotational  component  (w  x  r)  •  n,-  and  assume 
that  each  of  them  is  either  positive  or  negative,  there  are  2x2  =  4  combinations  of  signs. 
But  once  we  know  the  sign  of  the  flow  r  •  n,,  one  of  these  four  combinations  is  no  longer 
possible.  This  observation  has  been  used  in  the  development  of  global  constraints  for  3D 
motion  estimation.  Choosing  directions  n,  in  particular  ways  the  signs  of  (r  •  n,)  form  global 
patterns  of  positive  and  negative  areas  on  the  image  [3-5].  These  patterns,  whose  location 
and  form  encodes  information  about  3D  motion,  were  successfully  used  in  the  recovery  of 
egomotion.  In  this  paper,  by  pursuing  a  theoretical  investigation  of  the  amount  of  infor¬ 
mation  present  in  directional  flow  fields,  we  demonstrate  the  power  of  the  qualitative  image 
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measurements  already  used  empirically,  and  justify  their  utilization  in  global  constraints  for 
three-dimensional  dynamic  vision  problems. 

The  organization  of  this  paper  is  as  follows:  In  Section  3  we  develop  the  preliminaries 
constraints  that  will  be  used  in  the  uniqueness  analysis.  Given  two  rigid  motions,  we  study 
what  the  constraints  are  on  the  surfaces  in  view  for  the  two  motion  fields  to  have  the  same 
direction  at  every  point.  From  these  constraints,  we  investigate  for  which  points  of  the  image 
one  of  the  surfaces  must  have  negative  depth.  The  locations  where  negative  depth  occurs 
are  described  implicitly  in  the  form  of  constraints  on  the  signs  of  functions  depending  on  the 
image  coordinates  and  the  two  three-dimensional  motions.  The  existence  of  image  points 
whose  associated  depth  is  negative  ensures  that  the  two  rigid  motions  cannot  produce  motion 
fields  with  the  same  direction.  In  Section  4,  which  contains  the  main  uniqueness  proof,  we 
study  conditions  under  which  two  rigid  flow  fields  could  have  the  same  direction  at  every 
point  on  a  half  sphere  (i.e.,  conditions  under  which  there  do  not  exist  points  of  negative 
depth),  and  we  visualize  the  locations  of  negative  depth  on  the  sphere.  Section  5  is  devoted 
to  the  treatment  of  special  cases.  As  a  byproduct  of  the  analysis,  in  Section  6  we  investigate 
the  ambiguity  of  rigid  motion  for  full  flow  assuming  that  depth  has  to  be  positive,  and  show 
that  any  two  different  motions  can  be  distinguished  on  a  hemispherical  image  from  full  flow. 
Section  7  summarizes  the  results.  Appendix  A  studies  whether  more  than  two  rigid  motions 
could  produce  the  same  directional  flow  field,  and  the  rest  of  the  Appendices  (B-F)  describe 
and  prove  a  number  of  geometric  properties  used  in  the  main  part  of  the  paper. 

3  Critical  Surface  Constraints 

Let  us  assume  that  two  different  rigid  motions  yield  the  same  direction  of  flow  at  every  point 
in  the  image.  Let  ti  and  ui  be  translational  and  rotational  velocities  of  the  first  motion, 
and  let  and  u}2  be  translational  and  rotational  velocities  of  the  second  motion.  Since  from 
the  direction  of  flow  we  can  only  recover  the  directions  of  the  translation  and  rotation  axes, 
we  assume  all  four  vectors  ti,  t2,  and  a>2  to  be  of  unit  length.  Let  Zi{r)  and  Z2{t)  be  the 
functions,  mapping  points  r  on  the  image  into  the  real  numbers,  that  represent  the  depths 
of  the  surfaces  in  view  corresponding  to  the  two  motions.  In  the  future  we  will  refer  to  Zi 
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and  Zi  as  the  two  depth  maps.  In  this  section  we  investigate  the  constraints  that  must  be 
satisfied  by  Z\  and  Z^  in  order  for  the  two  fiow  fields  to  have  the  same  direction. 

We  assume  that  the  two  depths  are  positive,  and  allow  Z\  or  Z2  to  be  infinitely  large. 
Thus  we  assume  IjZx  >0  and  I/.Z2  >  0. 

3.1  Notation 

We  start  by  defining  some  notation: 

/^(r)  =  [wiCJsr] 

ft{r)  =  [titar]  (2) 

gij{Y)  =  (iVi  X  r)  ■  (tj  X  r)  for  i,  j  =  1,2 

where  [abc]  =  (a  x  b)  •  c  denotes  the  triple  product  of  vectors  a,  b  and  c. 

These  functions  have  a  simple  geometric  meaning.  If  Wi  x  a;2  =  0,  then  /a;(r)  =  0  for 
any  r.  If  u>i  x  u}2  ^  0,  then  fw{^)  is  zero  for  points  r  lying  on  a  geodesic  passing  through 
(jDi  and  ^2-  In  this  case  /o;(i')  =  0  defines  the  locus  of  points  r  where  Uroti(r),  the  rotational 
component  of  the  first  motion,  is  parallel  to  Urot2(j')5  the  rotational  component  of  the  second 
motion.  Similarly  ft{r)  is  either  zero  everywhere,  or  it  is  zero  for  points  lying  on  a  geodesic 
passing  through  ti  and  t2.  In  this  case  /t(r)  =  0  is  the  locus  of  points  r  where  Utri(i*)?  the 
translational  component  of  the  first  motion,  is  parallel  to  Utr2(r),  the  translational  component 
of  the  second  motion. 

If  uJi  =  0,  or  tj  —  0,  then  gij{r)  =  0  for  any  r.  If  they  are  non-zero,  then  gij{r)  is 
zero  at  points  lying  on  a  second  order  curve  consisting  of  two  closed  curves  on  the  sphere, 
the  so-called  zero  motion  contour  of  motion  (tj,a;i).  Equation  gij{Y)  =  0  defines  the  locus 
of  points  where  Urot.(r)  is  parallel  to  Utrj(r)  (see  Appendix  B).  Throughout  the  paper  the 
functions  /a;(r),  /t(r),  gij{r)  and  the  curves  defined  by  their  zero  crossings  will  play  very 
important  roles. 

To  simplify  the  notation  we  will  usually  drop  r  and  write  only  /,•  and  gij  where  the  index 
i  in  fi  can  take  values  t  and  u).  There  is  a  simple  relationship  between  fi  and  gij.  Let 

ui  =  (r  X  (ti  X  r))  X  (wi  x  r)  =  -[wi  r  (ti  x  r)]r  =  -gu  r 
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U2  =  (r  X  (t2  X  r))  X  (a;2  x  r)  =  -g22  r 
Since  we  assume  (r  •  r)  =  1,  we  obtain 

ui  •  U2  =  giig22  (3) 

Let 

U3  =  (t2  X  r)  X  (r  X  (a;2  x  r))  =  -U2 

Then 

Uj  •  U2  =  — Ui  •  U3  =  ftfuj  +  gi2g2i  (4) 

From  equations  (3)  and  (4)  we  get 

911922  =  ftfu  +  912921  (5) 


3.2  Conditions  for  ambiguity 

Assuming  that  motion  (ti,a;i)  with  depth  map  Zi  and  motion  (t2,a?2)  with  depth  map  Z2 
give  rise  to  flow  flelds  with  the  same  direction  at  any  point;  then  there  exists  g  >  0  such 
that 

-  4-(r  X  (ti  X  r))  -  wi  X  r  =  (-^(r  x  (ts  x  r))  -  wa  x  r)  (6) 

Zi  ^2 

By  projecting  the  vector  equation  (6)  on  directions  ta  x  r  and  r  x  (wa  x  r)  we  obtain  two 
scalar  equations 


^[titar]  +  (wi  X  r)  •  (ta  X  r)  =  /i  (wa  x  r)  •  (ta  x  r) 

Z\ 

■^(ui2  X  r)  •  (ti  X  r)  -  [wiwar]  =  g  -^{<^2  x  r)  •  (ta  x  r) 
Z\  m 

Since  g  is  positive,  from  (7)  and  (8)  we  get  constraints  on 

sgn(^/<  +5ri2)  =  sgn(5r22) 

sgn(^5f2i  -  U)  =  sgn(^522) 

Zi  Z2 

where  sgn(-)  denotes  the  sign  function. 


(7) 

(8) 

(9) 

(10) 
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Let  us  define  si  =  —g^lft  and  5^  =  fulg2\-  At  any  point,  /,■  and  gij  are  constant,  so 
equations  (9)  and  (10)  provide  simple  constraints  on  We  call  them  the  si-constraint  and 
the  Sj-constraint  respectively. 


Similarly  we  can  project  equation  (6)  on  vectors  ti  x  r  and  r  x  (uji  x  r)  and  obtain 


constraints  on 


sgn(^ii)  =  sgn{-^  ft  +  g2i) 

sgn(-^^ii)  =  sgn(^^i2  +  f^) 

Zj\  Zj2 


(11) 

(12) 


We  define  =  g2\l fu  ^2  =  —fwlgu-  Equations  (11)  and  (12)  provide  constraints  on 
and  we  thus  call  them  the  S2-constraint  and  the  s^-constraint. 


Let  us  now  interpret  these  constraints:  ^  is  always  non-negative;  thus,  if  the  two  motions 
(ti,cji),  (t2,tti2)  with  their  corresponding  depth  maps  Z\  and  Z2  produce  flow  with  the  same 
direction,  the  depth  Z\  must  satisfy 


either  ^Jt  +  5^12  >  0  and  ^g2i  -  fw  >  0 
or  ^ft  +  gi2<0  and  ^p2i  -  fw  <  0 

Thus  Zi  has  a  relationship  to  the  surfaces: 


and 


Z(r)  = 
Z(r)  = 


1 


5i(r) 

1 

5i(r) 


(13) 

(14) 


Equations  (13)  and  (14)  provide  hybrid  definitions  of  scene  surfaces.  To  express  the  surfaces 
in  scene  coordinates  R,  we  substitute  in  the  above  equations  Z{r)v  =  R.  Dividing  (13)  by 
Z{r)  and  replacing  Z{rY  by  R^  in  (14)  we  obtain 


(ti  X  t2)  •  R  +  (wi  X  R)  •  (t2  X  R)  =  0  (15) 

and  (012  X  R)  •  (ti  x  R)  —  ((aii  x  012)  •  R)R^  =  0  (16) 


Thus  we  see  that  Z\  is  constrained  through  (15)  by  a  second  order  surface  and  through  (16) 
by  a  third  order  surface.  At  some  points  it  has  to  be  inside  the  first  surface  and  at  some 
points  it  has  to  be  outside  the  first  surface.  In  addition,  at  some  points  it  has  to  be  inside  the 
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Figure  3:  Two  rigid  motions  (ti,a>i),  constrain  the  possible  depth  Zi  of  the  first 

surface  by  a  second  and  a  third  order  surface.  The  particular  surfaces  shown  in  the  coor¬ 
dinate  system  of  the  imaging  sphere,  projected  stereographically,  correspond  to  the  motion 
configuration  of  Figure  7. 

second  surface  and  at  some  points  it  has  to  be  outside  the  second  surface.  Figure  3  provides 
a  pictorial  description  of  the  two  surfaces  constraining  Zi. 

Analogous  to  the  above  derivation,  from  equations  (11)  and  (12)  we  obtain  a  further 
second  and  third  order  surface  pair,  which  constrain  the  depth  map  Z2- 

3.3  Interpretation  of  surface  constraints 

We  next  describe  the  si-  and  s^-constraints  in  detail.  For  convenience,  we  express  these 
constraints  for 

If  ft  =  0,  then  the  srconstraint  is  sgii{gi2)  =  sgn(5'22)?  i-e-?  it  does  not  depend  on  Zi. 
Thus  it  is  either  satisfied  by  any  Zi,  or  it  cannot  be  satisfied  by  any  Zi. 

If  ft  ^  0,  then  we  get  the  si-constraint  sgn(;|j-  —  si)  =  sgn(5'22)sgn(/t).  So  the  Si- 
constraint  is 

•  ;^  >  5i,  if  sgn(c?22)  sgn(/t)  >  0 

•  ^  =  si,  if  sgn(^r22)  =  0 

•  <si,  if  sgn(5'22)  sgn(/t)  <  0 

If  ^21  =  0,  then  the  s'^-constraint  does  not  depend  on  Zi. 


If  5^21  7^  0,  we  get  sgn(^  -  =  sgn(^5r22)sgn(£f2i).  Since  we  assume  ;^  >  0,  sgn(^^r22) 

is  either  0,  or  sgn(5r22)-  So  the  s'^ -constraint  is 

•  ^  >  if  sgn(^22)sgn(<72i)  >  0 

•  i  if  sgn(^22)  =  0 

•  if  sgn(<?22)sgn(^2i)  <  0 

At  each  point  we  have  the  additional  constraint  >  0.  If  all  three  constraints  can 
be  satisfied  simultaneously  at  a  point  in  the  image,  then  there  is  an  interval  (bounded  or 
unbounded)  of  values  of  Zi  satisfying  them.  If  the  constraints  cannot  be  satisfied,  this  means 
that  the  two  flows  at  this  point  cannot  have  the  same  direction  and  we  say  that  we  have  a 
contradictory  point. 

In  the  following  table  we  summarize  the  three  constraints  on  According  to  the 
inequality  relationships  from  the  si-constraint  and  the  Sj-constraint  we  classify  the  image 
points  into  four  categories  (type  I-IV).  The  table  analyzes  the  general  case,  at  a  point  where 
fi  7^  0  and  gjj  7^  0. _ 


Type 

sj-constraint 

Sj-constraint 

-f-  solution  interval 

Solution  exists  if 

I 

1/^1  ^ 

1/Zi  >  si 

(max(si,si,  0),  00) 

always 

II 

IjZi  >  Si 

1/^1  <  5i 

(max(si,0),si) 

si  >  0  and  Si  <  si 

III 

1/Zi  <  Si 

1/^1  >  5i 

(max(si,0),si) 

Si  >  0  and  si  <  Si 

IV 

IjZi  <  Si 

l/Zi  <  s[ 

(0,min(si,si)) 

Si  >  0  and  si  >  0 

If  some  of  fi,  pij  are  zero  at  a  point,  we  may  obtain  constraints  that  do  not  depend  on 
Z\,  or  equality  constraints. 

In  the  table  above  each  image  point  is  assigned  to  one  of  four  categories  (see  Figure  4  for 
an  example).  Whether,  for  a  given  image  point,  there  actually  exists  a  value  for  ^  satisfying 
the  constraints  depends  on  whether  the  solution  interval  at  that  point  is  empty  or  not.  Thus 
we  classify  all  image  points  on  the  sphere  into  three  categories.  A,  B,  and  C,  depending  on 
the  kind  of  solution  interval  that  exists  for  Z\:{K)  there  exists  no  solution  for  Zi,  or  there 
exists  a  solution,  and  the  interval  for  Zi  is  (B)  bounded,  or  (C)  unbounded.  In  the  latter 
two  cases,  we  can  also  check  whether  the  interval  has  a  lower  bound  greater  than  0. 
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Figure  4:  Classification  of  image  points  by  si-constraint  and  5'1-constraint. 


The  clcissification  of  a  point  into  one  of  the  categories  (I-IV)  depends  on  the  signs  of  /t, 
^21,  and  g22-  The  existence  of  a  solution  interval  at  a  point  also  depends  on  the  signs  of 
and  gx2  at  that  point  and  also  on  the  relative  values  of  5i  and  s'l,  i.e.,  on  the  sign  of  5i  —  s[. 

Functions  /,(r)  and  gijir)  are  polynomial  functions  of  r.  To  find  out  where  they  change 
sign,  it  is  enough  to  find  points  where  they  are  zero.  The  sign  of  5i  —  Sj  is  more  complicated, 
since  Si(r)  and  5'i(r)  do  not  have  to  be  continuous.  However  their  discontinuities  occur  at 
points  where  /t(r)  =  0  or  =  0.  Thus  sgn(5i  —  s^)  can  change  at  those  points  and  at 

points  where  Si  -  si  =  0.  Using  (5),  we  can  write 


5i  —  si 


gl2  fw 

ft  921 


1 

ft  921 


{ftfw  +  9l292\) 


9\1 922 
ft  921 


(17) 


So  we  see  that  sgn(si  —  si)  can  change  only  at  points  where  at  least  one  of  /i,  gij  is  zero. 
At  points  where  g22  —  0  have  the  si-constraint  ^-^d  the  si-constraint 

From  equation  (17)  we  obtain  si  =  si;  thus  at  these  points  the  depth  Zi  is  uniquely  defined. 

Let  us  consider  the  implicit  curves  /j(r)  =  0  and  9ij{r)  =  0.  In  the  general  case,  these 
equations  describe  two  geodesics  and  four  zero  motion  contours.  Each  of  the  curves  divides 
the  sphere  into  areas  where  the  solution  interval  for  Zi  could  be  different  (areas  of  class  A, 
B,  C).  However,  not  every  point  on  the  curves  separates  different  areas.  Inside  any  of  the 
areas,  all  the  points  have  the  same  classification  (for  example  an  infinite  solution  interval  for 
Zi  with  a  positive  lower  bound).  Figure  5  shows  an  example  of  this  classification,  although 
the  derivation  of  how  to  actually  obtain  the  areas  where  there  does  not  exist  a  solution  is 


deferred  to  the  next  subsection. 


Figure  5:  Classification  of  image  points  according  to  the  solution  interval.  (The  correspond¬ 
ing  motion  configuration  is  displayed  in  Figure  7.) 

Up  to  this  point  we  have  been  discussing  only  the  constraints  for  Z\.  Similarly,  from  (11) 
and  (12)  we  have  at  any  point  the  S2-constraint  (^  >  si-,  ^  —  52,  or  <  S2),  and  the 
S2-constraint  (^  >  s'^-,  ^  =  S2,  or  ^  <53).  We  obtain  the  same  curves  dividing  the  sphere 
into  areas  such  that  inside  any  of  the  areas,  the  type  of  solution  interval  for  Z2  is  the  same. 

Now  we  can  summarize  the  results.  The  curves  /i(r)  =  0  and  ^ij(r)  =  0  separate  the 
sphere  into  a  number  of  areas.  Each  of  the  areas  is  either  contradictory  (i.e.,  containing  only 
contradictory  points),  or  ambiguous  (i.e.,  containing  points  where  the  two  motion  vectors 
can  have  the  same  direction).  Two  different  rigid  motions  can  produce  ambiguous  directions 
of  flow  if  the  image  contains  only  points  from  ambiguous  areas.  There  are  also  two  scene 
surfaces  constraining  depth  Z\  and  two  surfaces  constraining  depth  Z2.  If  the  depths  do  not 
satisfy  the  constraints,  the  two  flows  are  not  ambiguous. 

3.4  Contradictory  points 

In  this  section  we  investigate  conditions  that  must  be  satisfied  when  a  point  is  contradictory. 
Since  the  type  of  solution  for  Zx  and  Z2  depends  on  the  signs  of  fi  and  gij,  we  want  to 
describe  sign  combinations  that  yield  a  contradiction.  We  investigate  the  general  case,  i.e., 
we  assume  fi  ^  0,  and  gij  7^  0  and  use  the  resulting  constraints  in  Section  4.  Special  cases 
are  treated  separately  in  Section  5. 
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There  are  two  simple  conditions  yielding  contradiction  for  Zi,  one  for  the  si-constreiint 
and  one  for  the  s'l-constraint.  There  is  no  solution  for  Zi  li  ^  <  5i  and  si  <  0.  This 
happens  under  the  following  condition  Ci. 

sgn(/t)  =  sgn(5^i2)  =  -sgn(5'22)  (18) 

which  is  derived  from  equation  (9).  Similarly,  from  (10)  we  get  a  contradiction  if  <  -si 
and  s'l  <  0,  i.e.,  under  condition  C2'- 

sgn(^)  =  -sgn(i?2i)  =  sgn(5f22)  (19) 

We  get  similar  conditions  for  Z2.  There  is  no  solution  for  Z2,  if  -^  <  S2  and  S2  <  0,  or  if 
^  <  S2  and  S2  <  0.  This  happens  under  conditions  C3  and  C4: 

sgn{  ft)  =  -sgn(5f2i)  =  sgn(5^ii)  (20) 

and 

sgn(/a,)  =  sgn(5^i2)  =  -sgn(5^ii)  (21) 

We  call  these  four  constraints  (Ci-C'4)  Contradictory  Point  conditions,  or  CP-conditions  for 
short.  Next  we  show  that  a  point  (where  /i  #  0  and  gtj  ^  0)  is  contradictory  if  and  only  if 
at  least  one  of  the  four  conditions  is  satisfied. 

Let  us  assume  that  conditions  (18)  and  (19)  are  not  satisfied  at  some  point,  but  we  have 
a  contradiction  for  Zi.  Then  the  point  must  be  of  type  II  or  III,  since  there  is  always  a 
solution  for  points  of  type  I,  and  a  point  of  type  IV  is  contradictory  only  if  (18)  or  (19) 
holds. 

For  a  point  of  type  II,  ^  <  s'j,  but  (19)  is  not  satisfied,  so  we  have  s'l  >  0.  A  contradiction 
is  possible  only  if  Si  >  This  happens  when  sgn(5f22)  =  sgn(/t),  sgn(5'22)  sgn(5f2i), 
sgn{gi2)  sgn(/t)  and  sgn(^)  =  sgn(5'2i),  i-e.,  when 

sgn(5f22)  =  sgn(/t)  =  -sgn(5'2i)  =  -sgn(5'2i)  =  -sgn(/a,) 

and  si  -  s[  >  0.  Since  sgn(/t)  =  -sgn(5'2i),  from  (17)  we  obtain 

sgn(5'n)  =  sgn(5f22)  =  sgn(/()  =  -sgn(5f2i) 
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Thus  in  this  case  condition  (20)  holds. 

We  obtain  the  same  result  for  points  of  type  III.  Since  (18)  is  not  satisfied,  we  have 
Si  >0.  So  a  contradiction  is  possible  only  if  Si  <  s^.  This  happens  when 

sgn(5'22)  =  -sgn(/i)  =  sgn(p2i)  =  sgn(^2i)  =  sgn(/^) 
and  Si  —  s'l  <  0.  Since  sgn(/t)  =  — sgn(5f2i),  we  obtain 

sgn(£'ii)  =  -sgn(p22)  =  sgn(/i)  =  -sgn(p2i) 
and  again  condition  (20)  holds. 

Thus  if  there  is  no  solution  for  Z\,  at  least  one  of  conditions  (18),  (19)  and  (20)  must 
hold.  Similarly  if  there  is  no  solution  for  Z2,  at  least  one  of  conditions  (18),  (20)  and  (21) 
must  hold. 

By  examination  of  all  the  possibilities,  we  can  show  that  at  any  point,  either  none  of 
the  CP-conditions  holds  (and  the  point  is  ambiguous),  or  exactly  two  of  the  conditions  hold 
(and  the  point  is  contradictory). 

3.5  Antipodal  pairs  of  points 

In  this  section  we  investigate  constraints  for  a  point  r  and  its  antipodal  point  — r  to  be  both 
ambiguous  or  to  be  both  contradictory. 

Again  we  describe  a  general  case,  i.e.,  assume  fi^O  and  gij  ^  0.  We  have  /t(— r)  = 
-/t(r),  /a;(-r)  =  -U{r),  but  gij{-r)  =  gij{T).  If  sgn(5ri2(r))  ^  sgn(5-22(r)),  then  con¬ 
dition  (18)  holds  either  at  r,  or  at  — r.  We  get  similar  results  for  the  remaining  three 
CP-conditions.  Thus  both  point  r  and  point  — r  are  ambiguous  only  if 

sgn(pii(r))  =  sgn(5ri2(r))  =  sgn(5r2i(r))  =  sgn(5-22(r))  (22) 

Point  r  and  point  — r  can  also  both  be  contradictory.  As  shown  in  Appendix  F,  this 
happens  when 

sgn(pii(r))  =  sgn(p22(r))  =  -sgn(p2i(r))  =  -sgn(^22(r))  (23) 
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4  The  Geometry  of  the  Depth-Positivity  Constraint 

In  the  last  section  we  found  that  if  the  CP-conditions  hold  at  a  point  on  the  imaging  surface, 
then  one  of  the  depth  values  has  to  be  negative  and  thus  the  point  is  contradictory.  In  this 
section  we  investigate  these  constraints  further;  in  particular  we  would  like  to  know  under 
what  conditions  two  rigid  motions  cannot  be  distinguished  if  our  imaging  surface  is  a  half 
sphere  or  an  image  plane,  and  we  are  interested  in  studying  and  visualizing  the  locations  of 
areas  where  the  CP-conditions  are  met. 

Considering  as  imaging  surface  the  whole  sphere,  two  different  rigid  motions  cannot 
produce  flow  of  the  same  direction  everywhere.  As  shown  in  Section  3.5,  two  antipodal 
points  r  and  — r  are  ambiguous  only  if  (22)  holds.  Thus  for  any  point  on  curve  gij  =  0,  since 
the  sign  of  gij  is  positive  on  one  side  of  the  curve  and  negative  on  the  other,  there  must  exist 
a  neighborhood  either  around  r  or  around  -r  where  there  is  a  contradiction. 

We  are  now  ready,  using  the  machinery  already  developed,  to  study  uniqueness  properties. 
As  in  the  previous  section,  we  assume  that  vectors  ti,  t2,  Wi,  and  UI2  are  of  unit  length. 

4.1  Half  sphere  image:  The  general  case 

Let  us  assume  that  the  image  is  a  half  of  the  sphere.  Let  us  also  assume  that 

(wi  X  W2)  •  (ti  X  t2)  0  (24) 

We  show  that  under  this  condition  the  two  rigid  motions  cannot  produce  motion  fields  with 
the  same  direction  everywhere  in  the  image. 

Let  us  consider  the  projections  of  and  on  a  geodesic  n  connecting  ti  and  t2. 
Projection  onto  the  geodesic  is  well  defined  for  all  points  r  such  that  r  x  (ti  x  12)  0.  Since 

we  assume  (24),  the  projections  of  both  u>i  and  u>2  are  well  defined.  The  proof  is  given  in 

parts  A  and  B. 

A:  Let  us  first  assume  that  one  of  Wi,  a;2  does  not  lie  on  geodesic  n.  Without  loss  of 
generality,  let  it  be  wi,  i.e.,  let  [tit2Wi]  7^  0. 
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Figure  6:  Possible  sign  combinations  of  gn  and  in  the  neighborhood  of  ri. 

The  projection  of  u)i  onto  n  is 

ri  =  ±(ti  X  t2)  X  {oj-i,  X  (ti  X  t2))  (25) 

where  the  sign  is  chosen  so  that  ri  is  in  the  image.  Then 

/^(ri)  =  ±((wi  X  (J}2)  X  (ti  X  t2))  •  (wi  X  (ti  X  t2)) 

(26) 

—  T[ti  t2  u;i]((ct;i  x  ^2)  •  (ti  x  t2))  ^  0 

and  5fii(ri)  =  5'i2(ri)  =  0.  So  the  s^-constraint  is  sgn(;^0)  =  sgn(;^0  +  /„)  or  =  0. 
Clearly,  this  constraint  cannot  be  satisfied,  so  ri  is  a  contradictory  point. 

We  can  also  show  that  at  least  one  of  the  areas  around  point  ri  is  contradictory.  Point  ri 
lies  on  zero  motion  contours  5rii(r)  =  0  and  i?i2(r)  =  0.  If  the  two  contours  cross  at  this  point 
(Appendix  C  shows  that  5rii(r)  =  0  and  5'i2(r)  =  0  at  ri  cannot  be  tangent),  we  obtain  four 
areas  in  the  neighborhood  of  ri,  and  all  four  possible  sign  combinations  of  and  gi2-  If  we 
look  at  points  close  enough  to  ri  (so  that  does  not  change  sign),  then  condition  (21)  is 
satisfied  in  one  of  the  areas,  and  that  area  is  contradictory.  For  an  illustration  see  Figure  6. 

B:  Now  we  need  to  consider  the  situation  where  both  and  0^2  lie  on  geodesic  n,  i.e., 
[tit2U>i]  =  [tit2  W2]  =  0.  Let  us  consider  point  u>i.  We  know  /t(a;i)  =  gii{uJi)  =  ^'12(^1)  = 
0.  Also  g2i{(*^i)  —  {<-^2  X  Wj)  •  (tj-  X  Wi).  Since  both  coi  and  UJ2  lie  on  geodesic  n,  (a;2  x  oJi) 
is  parallel  to  (t,-  x  Wi).  Thus  g2i{<’^i)  is  zero  only  if  (tj-  x  Wi)  is  zero.  However,  from  (24)  we 
have  (ti  X  t2)  /  0,  so  either  g2i  or  g22  is  non- zero  at  Wi. 
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If  0,  then  condition  (11)  cannot  be  satisfied  and  oji  is  a  contradictory  point. 

Again,  it  is  not  a  singular  point.  The  line  tangent  to  at  has  direction  Wi  x  ti  (and 
ct?!  X  ti  7^  0,  since  ^21(^1)  /  0)?  so  pn  is  perpendicular  to  n  at  this  point.  Since  /i  is 
identical  to  n,  curves  gw  =  0  and  ft  =  0  create  four  areas  around  uji  with  all  possible  sign 
combinations.  Thus  in  one  of  the  areas,  condition  (20)  holds,  and  we  obtain  a  contradictory 
area. 

If  ^'22(^1)  7^  0,  then  condition  (9)  cannot  be  satisfied  at  u>i.  Again,  at  least  one  area 
around  oJi  is  contradictory,  since  contour  gu  =  0  is  perpendicular  to  n  at  this  point.  This 
concludes  the  proof  that  if  (24)  is  satisfied  there  exist  contradictory  areas  on  the  half  sphere. 
Section  4.2  discusses  the  case  when  (24)  is  not  satisfied. 

The  rest  of  this  section  describes  properties  of  the  contradictory  areas  in  order  to  provide 
a  geometric  intuition. 

Just  as  we  projected  Wi  on  geodesic  n  connecting  ti  and  t2  to  obtain  ri,  we  project  UJ2 
on  n  to  obtain  r2,  and  we  project  tj  and  t2  on  geodesic  /,  connecting  uJi  and  0^2,  to  obtain 
r3  and  r4  (see  Figure  7).  Point  r2  is  at  the  intersection  of  ft  =  0,  P21  =  0,  and  g22  =  0;  rs  is 
at  the  intersection  of  /u,  =  0,  5-11  =  0,  and  521  =0;  and  r4  is  at  the  intersection  of  /u,  =  0, 
gi2  =  0,  and  P22  =  0-  By  the  same  argumentation  as  before,  at  each  of  the  points  we  can 
choose  two  of  the  contours  fi  =  0  and  gij  =  0  passing  through  the  point  and  we  obtain  four 
areas  of  different  sign  combinations  in  the  corresponding  terms  fi  and  gij  around  the  point; 
it  can  be  shown  that  one  of  these  areas  is  contradictory  because  one  of  the  CP-conditions  is 
met. 

The  CP-conditions  are  constraints  on  the  signs  of  the  terms  fi  and  gij.  Thus  the  bound¬ 
aries  of  the  contradictory  areas  are  formed  by  the  curves  fi  =  0  and  gij  =  0.  As  we  have 
shown  the  contradictory  area  and  its  boundaries  must  contain  the  points  ri,  r2,  ra,  and  r4. 
For  some  motion  configurations  the  boundaries  also  might  contain  ti,  t2,  Wi,  and  u>2-  If 
can,  however,  be  verified  that  no  neighborhood  around  ti,  t2,  Wi,  and  W2  needs  to  be  con¬ 
tradictory.  It  can  also  be  verified,  by  examining  all  the  possibilities  for  the  signs  of  terms  /, 
and  gij  in  the  CP-conditions,  that  points  ti,  t2,  and  UP2  cannot  lie  inside  a  contradictory 
area,  since  at  least  one  of  their  neighboring  areas  is  ambiguous.  Figures  8  and  9  show  the 
contradictory  areas  for  both  halves  of  the  sphere  for  two  different  motion  configurations. 
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Figure  7:  Separation  of  the  sphere  through  curves  fi  =  0  and  gij  =  0.  Each  of  tj,  t2,  Wi, 
^2,  Ti,  1^2)  1*3,  and  r4  lies  at  the  intersection  of  three  curves. 


(a)  (b) 

Figure  8;  Contradictory  areas  for  both  halves  of  the  sphere  for  the  two  motions  shown  in 
Figure  7. 

Finally,  let  us  consider  the  boundaries  of  the  contradictory  areas.  As  defined  in  Section  3, 
we  allow  the  depths  of  the  surfaces  in  view  to  take  any  value  greater  than  zero  (including 
infinity).  Thus  at  any  point  r  the  motion  vector  f  could  be  in  the  direction  of  Urot(r),  but  not 
in  the  direction  of  ftr(r).  This  allows  us  to  describe  the  depth  values  at  possible  boundaries 
of  a  contradictory  area:  At  points  on  curve  fu,  =  0  both  Zi  and  Z2  can  be  infinite,  thus 
boundary  points  on  this  curve  are  not  elements  of  the  contradictory  area.  Boundary  points 
on  all  other  curves  (/<  =  0,  or  g^j  =  0)  are  contradictory,  since  one  of  the  depths  Zi  and  Z2 
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Figure  9:  (a)  Motion  configuration,  (b)  and  (c)  Contradictory  areas  for  both  halves  of  the 
sphere. 

would  have  to  be  zero. 

4.2  Half  sphere  image:  The  case  when  (ti  x  12)  is  perpendicular  to  (a?i  x  0^2) 

In  this  section  it  is  shown  that  there  could  exist  (ti,Wi)  and  (t2,W2),  with  (ti  x  t2)  perpen¬ 
dicular  to  (wi  X  a;2),  such  that  there  exist  no  contradictory  areas  in  one  hemisphere. 

First  we  investigate  possible  positions  of  points  tj,  t2,  and  W2  on  the  hemisphere, 
bounded  by  equator  q.  Then  we  describe  additional  conditions  on  the  orientation  of  vectors 
tj-  and  (jJi  with  respect  to  the  hemisphere. 

As  shown  in  Section  3.5,  two  antipodal  points  r  and  — r  can  be  ambiguous  only  if  (22) 
holds.  Thus  if  the  border  of  the  area  defined  by  (22)  intersects  q,  there  will  be  a  contradiction 
in  the  image. 

If  curve  gij  =  0  intersects  q  at  point  p,  at  least  one  of  the  areas  around  p  does  not  satisfy 
condition  (22).  Unless  ti,  t2,  and  UJ2  all  are  on  the  boundary  of  the  hemisphere  (and  then 
the  motions  are  not  ambiguous),  there  is  a  contradictory  area  in  the  image  (either  around 
p,  or  around  — p). 

Let  no  be  the  normal  to  the  plane  of  q.  By  intersecting  the  zero  motion  contour  gij  —  0 
with  the  border  q  of  the  hemisphere  (see  Appendix  D),  we  find  that  real  solutions  for  the 
intersection  point  are  obtained  only  if 

/  =  [twno]^  -  4(w  -  no)(t  •  no)(w  •  t)  >  0.  (27) 
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A  half  sphere  contains  for  each  of  the  translation  vectors  and  the  rotation  vectors 
uji,  exactly  one  of  the  vectors  +ti  or  — tj-  and  +a;j-  or  — Wj-.  Let  us  refer  to  the  vectors 
in  the  considered  hemisphere  as  ii  and  tDj-.  From  equation  (27),  taking  into  account  that 

(tj-  •  no)  >  0  and  (u;^  •  no)  >  0,  we  see  that  /  >  0,  either  if  for  any  tj,  (tDj-  •  tj)  <  0 

(i.e.,  Wi  and  tj  form  an  angle  greater  than  90°),  or  (w,  •  tj)  >  0  and  uji  and  tj  are  such 
that  [tju;,no]^  >  4(a;j  •  no)(tj  •  no)(a;i  •  tj),  which  means  that  w;  and  tj  must  be  close  to  the 
border. 

When  /i  =  0  is  perpendicular  to  =  0,  the  projections  of  wi  and  u}2  on  ft  =  0  and  the 
projections  of  ti  and  t2  on  =  0  coincide  in  one  point  ri,  i.e.,  ri  =  ta  =  rs  =  r4.  Point  rj 
lies  at  the  intersection  of  all  six  curves  fi  =  0  and  gij  —  0. 

Any  three  curves  fi  =  0,  gjj  =  0  and  g^i  —  0  (with  k  ^l)  intersect  only  in  rj  and  one  of 

the  points  ti,  t2,  ^i,  or  W2.  Furthermore,  since  all  the  zero  motion  contours  have  to  be  closed 
curves  on  the  hemisphere,  we  conclude  that  if  there  exists  a  contradictory  area,  it  also  has  to 
be  in  a  neighborhood  of  ri.  It  thus  suffices  to  consider  all  possible  sign  combinations  of  terms 
fi  and  gij  around  ri.  It  can  be  verified  that,  for  a  hemisphere  to  contain  only  ambiguous 
areas,  the  two  translations  have  to  have  the  same  sign,  that  is  sgn(ti  •  no)  =  sgn(t2  •  no).  Also 
the  two  rotations  have  to  have  the  same  sign,  i.e.,  sgn{u>i  ■  no)  =  sgn(w2  •  no).  Furthermore, 
the  relative  positions  of  ti,  t2,  Wi,  and  uj2  have  to  be  such  that 

sgn(((wi  X  u>2)  X  (ti  X  t2))  •  no)  =  sgn(ti  •  no)  sgn(wi  •  no) 

Intuitively  this  means,  when  rotating  =  0  in  the  orientation  given  by  the  rotations  in 
order  to  make  /^  =  0  and  ft  =  0  parallel,  then  the  order  of  points  ti  and  t2  on  /*  =  0  is 
opposite  to  the  order  of  points  Cji  and  0^2  on  f^,  =  Q  (moving  along  the  same  direction  along 
fu;  =  0  and  ft  =  0),  if  sgn(ti  •  no)  =  1.  Otherwise,  if  sgn(ti  •  no)  =  —1,  the  order  of  points 
— ti  and  — 12  on  ft  =  0  must  be  the  same  as  the  order  of  points  u>i  and  0^2  on  f^  =  0. 

In  summary,  we  have  shown  that  two  rigid  motions  could  be  ambiguous  on  one  hemi¬ 
sphere,  if  (ti  X  t2)  is  perpendicular  to  (o^i  x  0^2)5  but  only  if  certain  sign  and  certain  distance 
conditions  on  ti,  t2,  and  0^2  are  met.  In  addition,  as  shown  in  Section  3,  the  two  surfaces 
in  view  are  constrained  by  a  second  and  a  third  order  surface  (as  shown  in  equations  (15) 
and  (16)).  Figure  10  gives  an  example  of  such  a  configuration. 
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Figure  10:  Both  halves  of  the  sphere  showing  two  rigid  motions  for  which  there  do  not  exist 
contradictory  areas  in  one  hemisphere,  (a)  Hemisphere  containing  only  ambiguous  areas, 
(b)  Contradictory  areas  for  the  other  hemisphere. 

In  the  next  section  we  discuss  the  special  cases  and  show  that  they  do  not  allow  for 
ambiguity.  Thus  the  case  of  (tj  x  t2)  being  perpendicular  to  (ttJi  x  UJ2)  is  the  only  case  where 
two  motions  can  produce  the  same  direction  of  the  motion  field  on  a  hemisphere.  An  analysis 
concerned  with  ambiguities  due  to  more  than  two  rigid  motions  is  given  in  Appendix  A. 


5  Special  Cases 

In  previous  sections,  we  assumed  that  tj  x  t2  7^  0  and  ufi  x  W2  ^  0.  Here  we  show  that  if 
these  conditions  do  not  hold,  then  the  two  motions  are  not  ambiguous. 

In  Section  3  we  assumed  all  four  vectors  ti,  t2,  Wi,  and  0^2  to  be  of  unit  length.  Here  the 
four  vectors  can  also  be  zero.  Thus  we  have  two  different  motions  (i.e.  ti  ^  t2,  or  iVi  ^  072), 
such  that  t:  X  t2  =  0  or/and  u>i  x  0)2  =  0. 

To  cover  all  possible  cases  we  are  required  to  make  a  minor  assumption  about  the  depth 
Zi  and  Z2  for  the  case  where  (jOx  =  uj 2. 

Let  u>i  =  a;2-  Then  we  have  ti  ^  t2,  and  /u,  =  0  everywhere.  From  (10)  we  obtain  the 
constraint  sgn(^5f2i)  So  at  points  where  g2i  and  ^22  have  different  signs,  the 

only  possible  solution  is  =  0-  Infinite  values  for  both  depths  in  these  areas  would 
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result  in  pure  rotational  flow  fields  in  these  areas  and  thus  in  an  ambiguity.  The  same  kind 
of  ambiguity  would  occur  if  we  considered  the  full  flow.  Therefore  it  seems  reasonable  to 
assume  that  at  least  at  one  point  in  the  areas  where  gug22  <  0,  depths  Zi  and  Z2  are  not 
both  infinite.  Under  this  assumption  there  does  not  exist  ambiguity  for  the  case  of  (x>i  =  a?2. 
In  the  following  we  thus  assume  u>i  ^  a?2. 

Next  we  provide  a  lemma  that  will  be  of  use  in  the  following  proofs  concerned  with  special 
cases  as  well  as  in  the  proof  for  full  flow  in  Section  6. 


Lemma:  Let  uii  ^  0,  tj  ^  0.  As  in  the  previous  section,  let  the  image  be  a  half  sphere 
with  equator  g,  let  no  be  a  unit  vector  normal  to  the  plane  of  q.  Then  equation 


sgn(^/0  =  (Tsgn(^y) 


where  <7  =  ±1,  can  be  satisfied  everywhere  in  the  image  only  if  tj  x  no  =  0  and  cJj-  •  tj  =  0. 


Proof:  Since  and  tj  are  non-zero,  there  are  points  in  the  image  where  ^,j  ^  0.  Thus 
ti  X  t2  must  be  non-zero  and  geodesic  n  connecting  ti  and  t2  is  well  defined.  Equation  (28) 
can  be  satisfied  only  if  the  zero  motion  contour  is  degenerate,  i.e.,  w^-tj  =  0  (as  in  Figure  14b). 
Then  the  contour  consists  of  two  great  circles.  One  of  the  circles  must  be  identical  to  the 
geodesic  n,  and  the  other  circle  must  be  identical  to  §,  the  border  of  the  image.  This  is 
possible  only  if  tj  x  no  =  0  (see  Figure  11).  □ 


Figure  11:  If  sgn(^/t)  =  cr  sgn(5rij)  everywhere,  the  zero  motion  contour  gij  =  0  consists  of 
two  great  circles,  one  identical  to  the  border  of  the  hemisphere,  the  other  identical  to  ft  =  0. 

We  now  consider  two  special  cases  in  parts  A  and  B. 
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A;  Let  us  assume  that  all  t,-  and  w,-  are  non-zero. 

If  ti  X  t2  =  0,  then  ft  is  zero  everywhere.  Thus  from  condition  (9)  we  obtain  sgn(5'i2)  = 
sga{g22)-  Since  all  four  vectors  are  non-zero,  this  is  possible  only  if  Wj  x  0^2  =  0. 

So  we  only  need  to  consider  the  case  oJi  x  a;2  =  0,  and  since  we  also  assume  ^  u:2,  we 
have  0)2  =  —  wi.  Then  at  any  point  in  the  image,  g2j{f^)  —  ~gij{^)-  Thus  (9)  can  be  satisfied 
only  if  sgn(^/4)  =  sgn(5'22)-  According  to  the  lemma,  this  is  possible  only  if  t2  x  no  =  0. 

Similarly  (11)  can  be  satisfied  only  if  sgn{^  ft)  =  sgn(5r2i).  So  from  the  lemma  we  obtain 
ti  X  no  =  0.  Therefore  we  have  ti  x  t2  =  0,  function  ft  is  zero  everywhere,  and  the  motions 
are  contradictory. 

B:  If  one  of  the  motion  parameters  is  zero,  we  obtain  either  a  pure  translational  or  a  pure 
rotational  flow  field.  By  considering  all  the  possible  cases,  it  can  be  verified  that  the  two 
motions  are  not  ambiguous.  Here  we  just  consider  one  of  the  more  difficult  cases. 

Let  u>i  =  0,  u>2  ^  0,  ti  0,  and  t2  0.  Then  at  any  point,  gn  =  gi2  =  fw  =  0. 
So  from  (9)  we  obtain  sgn(;^/t)  =  sgn(5'22)5  from  (11)  we  have  0  =  sgn{—-^ft  +  gii),  ox 
sgn(^/t)  =  sgn(5'2i).  From  the  lemma,  this  is  possible  only  if  t2  x  no  =  0  and  ti  x  no  =  0; 
thus  again  we  obtain  ti  x  t2  =  0  and  the  motions  are  contradictory. 

If  two  of  the  motions  are  zero,  that  is  if  either  ti  =  t2  =  0  or  ti  =  W2  =  0  (or  equivalently 
t2  =  Wi  =  0)  we  obtain  either  two  rotational,  or  one  translational  and  one  rotational  field, 
which  obviously  cannot  have  the  same  direction. 

6  Ambiguities  of  the  Full  Flow 

Next  we  investigate  the  question  whether  there  can  be  any  ambiguities  at  all  if  we  consider 
the  complete  flow.  Horn  has  shown  in  [6]  that  two  motions  can  produce  ambiguous  flow 
fields  only  if  the  observed  surfaces  are  certain  hyperboloids  of  one  sheet.  We  show  that  if  we 
also  consider  the  depth  positivity  constraint  and  if  the  image  is  a  half  of  the  sphere,  then 
any  two  different  motions  can  be  distinguished. 

Let  the  image  be  a  hemisphere  bounded  by  equator  q.  Let  uq  be  a  unit  vector  normal  to 
the  plane  of  q.  As  in  [6],  let  us  assume  that  a  motion  (ti,Wi)  along  with  a  depth  map  Zi, 
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and  a  motion  (t2,W2)  along  with  a  depth  map  Z2,  yield  the  same  flow  field.  At  each  point 
we  obtain  a  vector  equation 

-  ^(r  X  (ti  X  r))  -  X  r  =  -i^(r  x  (t2  x  r))  -  U72  x  r  (29) 

Zj\  Zj2 

Projecting  on  directions  ti  x  r  and  t2  x  r,  we  obtain  equations  for  the  two  critical  surfaces 

•;^[tit2r]  +  X  r)  •  (t2  X  r)  =  0  (30) 

■^1 

^[tit2r]  +  (^a;  x  r)  •  (ti  x  r)  ==  0  (31) 

Z2 

where  =  W2  —  <^1- 

If  ti  ^  0,  t2  7^  0,  and  7^  0,  then  according  to  the  lemma  in  the  previous  section,  these 
equations  can  be  satisfied  everywhere  in  the  image  only  if  do?  •  ti  =  0,  tj  x  no  =  0,  6a;  =  0, 

and  t2  X  no  =  0.  Thus  we  we  obtain  ti  x  t2  =  0.  (This  case  corresponds  to  Section  4.5  in  [6], 
that  is,  to  the  case  when  both  critical  surfaces  consist  of  intersecting  planes.)  Therefore  we 
are  left  only  with  special  cases:  Ambiguity  can  occur  only  if  ti  x  t2  =  0,  or  Su  =  0. 

1^  ti  X  t2  =  0  and  6a;  7^  0,  from  constraint  (30)  we  get  for  any  r 

(6a;  X  r)  •  (t2  X  r)  =  0  (32) 

Since  6a;  7^  0,  t2  must  be  zero.  Similarly  from  constraint  (31)  we  get  ti  =  0.  Thus  we 
have  a  pair  of  rigid  motions  with  different  rotations  and  zero  translations.  Clearly  these  two 
motions  are  not  ambiguous. 

There  is  one  special  case  left,  6a;  =  0.  At  each  point  we  get  a  vector  equation 

-  ^(r  X  (ti  X  r))  =  -- ^(r  x  (ts  x  r))  (33) 

Zi  Z<2 

Since  we  have  two  different  motions  and  6a;  =  0,  we  know  ti  7^  t2.  So  the  equation  can  be 
satisfied  only  when  ^  ^  =  0  for  all  points  not  lying  on  geodesic  n  passing  through  ti 

and  t2.  If  we  do  not  allow  infinite  depth,  the  motions  are  not  ambiguous. 

7  Conclusions 

In  this  paper  we  have  analyzed  the  amount  of  information  inherent  in  the  directions  of  rigid 
flow  fields.  We  have  shown  that  in  almost  all  cases  there  is  enough  information  to  determine 

24 


up  to  a  multiplicative  constant  both  the  3D-rotational  and  3D-translational  motion  from 
a  hemispherical  image.  Ambiguities  can  result  only  if  the  surfaces  in  view  satisfy  certain 
inequality  and  equality  constraints.  Furthermore,  for  two  3D  motions  to  be  compatible  the 
two  translation  vectors  must  lie  on  a  geodesic  perpendicular  to  the  geodesic  through  the  two 
rotation  vectors.  With  this  analysis  we  have  also  shown  that  visual  motion  analysis  does  not 
necessarily  require  the  intermediate  computation  of  optical  flow  or  exact  correspondence. 
Instead,  many  dynamic  vision  problems  might  be  solved  with  the  use  of  more  qualitative 
flow  estimates  if  appropriate  global  constraints  are  found. 

Appendix 

Appendix  A  Ambiguity  due  to  more  than  two  motions 

In  this  appendix  we  investigate  whether  three  or  more  different  rigid  motions  and  their 
corresponding  surfaces  could  possibly  produce  the  same  direction  of  the  motion  field  on  a 
hemisphere.  We  present  proofs  contradicting  the  ambiguity  of  almost  all  combinations  of 
three  rigid  motions. 

Let  us  consider  any  three  different  rigid  motions  (ti,u;i),  (t2,ct’2)?  and  (t3,a>3),  such  that 
any  two  of  the  directional  motion  fields  produced  are  the  same,  i.e.  (t,  x  tj)  ■  {uJi  x  u^j)  =  0 
for  i  =  1, . . . ,  3,  i  =  1, . . . ,  3  with  i  ^  j.  In  the  following  proofs  it  will  be  shown  that  in 
general  there  exist  areas  in  the  image  where  the  corresponding  depth  Z3  cannot  at  the  same 
time  allow  motions  (ti,a>i)  and  (t3,a;3)  and  motions  and  (t3,a73)  to  produce  the 

same  directional  flow. 

Let  us  consider  the  intersections  of  the  zero  motion  contours  gu  =  0.  In  the  sequel  we 
consider  separately  in  part  A  the  general  case  where  two  of  the  zero  motion  contours  intersect 
in  at  least  two  points,  and  in  part  B  the  case  where  any  two  zero  motion  contours  are  tangent 
to  each  other.  (Appendix  E  describes  the  conditions  on  the  motion  parameters  for  two  zero 
motion  contours  to  be  tangential.) 

A:  Let  us  assume  that  two  of  the  zero  motion  contours  are  not  tangential;  let  these  be 
^11  =  0  and  g22  =  0.  Let  ri2  be  the  intersection  point  where  (ti  x  t2)  •  ri2  =  0  and 
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(u?i  X  u;2)  ■  ri2  =  0,  and  thus  Utri(ri2),  Vtr2(ri2),  ?^roti(ri2),  yrot2(i'i2)  are  parallel.  Let  pi2 
be  another  intersection  point  where  ffn  and  ^22  cross.  Vectors  Uroti(pi2)  and  Urot2(Pi2)  are 
not  parallel,  and  Utri(pi2)  =  — ^i^^roti  (P12)  and  Utr2(Pi2)  =  —  A2t^rot2(Pi2)  for  some  positive  Ai 
and  A2  (if  Ai  or  A2  were  negative,  point  pi2  would  be  contradictory).  Unless  5'33(pi2)  =  0, 
we  have  Utr3{pi2)  x  'Urot3(Pi2)  ^  0.  Figure  12a  shows  a  possible  configuration  of  the  motion 
vectors  at  pi2. 

We  next  consider  the  directions  of  Urot.  (i“)  and  Utri(r)  for  points  r  in  the  neighborhood 
of  pi2-  Let  no  be  a  unit  vector  in  the  direction  Utr3(pii)  x  Urot3(Pi2)-  For  i  =  1..2  the 
sign  of  (utr,(r)  x  ■yrot,(r))  •  no  changes  from  inside  ga  =  0  to  outside  ga  =  0  (that  is,  for 
example,  the  angle  between  Utri(r)  and  Uroti(r)  is  greater  than  180°  inside  gu  =  0  and  is 
smaller  than  180°  outside  gn  =  0,  or  vice  versa).  The  sign  of  (utr3(r)  x  t;i.ot3(r))  •  no  is 
the  same  in  a  sufficiently  small  neighborhood  around  pi2.  Since  gn  =  0  and  ^22  =  0  cross 
at  pi2  there  are  four  neighborhoods  around  pi2  with  all  four  possible  sign  combinations  of 
(utn  X  Uroti )  •  no  and  (utr2  x  Urotj)  •  no-  Thus  for  points  r  in  one  of  the  neighborhoods,  in  order 
for  Utn  (r)  +  Urotj  (r)  to  have  the  same  direction  as  Utr3(r)  +  Urot3(i'))  Z3  must  lie  in  an  interval 
[a,  6],  and  for  Utr2(r)  +  t'rot2(i‘)  to  have  the  same  direction  as  Utr3(r)  +  Vrot3(r),  Z3  must  lie 
in  an  interval  [c,  d],  where  the  intersection  of  [a,  6]  and  [c,  d]  is  empty.  Therefore  the  three 
motions  cannot  give  rise  to  the  same  direction  at  r.  For  an  example  see  Figure  12b. 

B:  We  next  consider  the  case  where  all  three  zero  motion  contours  are  tangent  to  each 
other.  For  the  case  where  not  all  three  are  tangential  at  the  same  point,  using  arguments 
similar  to  those  used  before,  we  prove  that  there  cannot  be  an  ambiguity. 

For  at  least  two  of  the  zero  motion  contours,  say  gn  =  0  and  522  =  0,  we  have  that 
at  the  intersection  point  ri2  the  two  curvature  vectors  «5ii(^i2)  of  =  0  2,nd  /C522(^i2) 
of  922  =  0  have  opposite  sign.  Also,  the  translational  and  rotational  components  are  such 
that  vtTi  =  AiUrot2  =  — ^2^’tr2  =  — ^s^roti  for  Ai,A2,A3  >  0  (see  Figure  13a  and  b  for  an 
illustration).  Let  no  be  a  unit  vector  in  the  direction  Utr3(ri2)  x  Urot3(ri2)-  At  points  r  in 
the  neighborhood  of  ri2  we  obtain  three  of  the  four  possible  sign  combinations  for  the  signs 
of  (t^roti(r)  X  t;tri(r))  •  no  and  (utr2(i’)  x  Urot2('f))  '  ^o-  In  both  areas  outside  gn  —  0  and 
outside  g22  =  0  we  have  ((t;roti(r)  x  Utri(r))  •  no)(utr2(r)  x  Ui.ot2(r))  •  no)  <  0,  but  in  one  of 
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(a)  (b) 


Figure  12:  (a)  Possible  motion  configuration  at  point  pi2.  (b)  There  must  exist  a  neighbor¬ 
hood  around  pi2  with  points  r,  such  that  for  no  =  Utr3(pi2)  X’^rot3(Pi2)/btr3(Pi2)  xUrot3(Pi2)|5 
(Uroti  X  Utri)  •  no  >  0  and  (utrj  X  Urotz)  •  no  <  0.  In  order  for  U3  =  utri(r)  +  Uroti(r)  = 
Ai(utr3(r)  +  Urot3(r)),  U3  has  to  be  in  the  sector  Si  and  Z3  has  to  take  values  in  the  interval 
(0,  b].  In  order  for  U3  =  t;tr2(r)  +  =  A2(utr3(r)  +  UrotjCr)),  U3  has  to  be  in  the  sector 

S2  and  Z3  has  to  take  values  in  the  interval  [c,  00]  with  b  <  c. 

the  areas  (utr2(r)  x  'yroti(i'))  •  no  >  0  and  in  the  other  (yirzCr)  x  Uroti(r))  •  no  <  0.  Since 
(utr3(r)  X  UrotsCr))  -  no  doesn’t  change  sign  in  the  neighborhood  of  ri2,  in  one  of  the  two  areas 
the  depth  Z3  of  the  third  surface  cannot  be  compatible  with  both  the  first  and  the  second 
motion  (see  Figure  13c). 


(a)  (b)  (c) 


Figure  13:  (a)  Intersection  of  zero  motion  contours  gu  =  0  and  ^22  =  0  at  point  ri2  with 
<  0-  (b)  Possible  motion  configuration  at  point  ri2.  (c)  At  point  r  in  one  of 
the  areas  outside  =  0  and  outside  5^22  —  0  (here  S2)  the  depth  of  Z3  cannot  be  compatible 
with  both  (utri(r),Uroti(r))  and  (utr2(r),Urot2(r))- 

Thus,  in  summary  we  have  shown  that  more  than  two  different  rigid  motions  can  hardly 
ever  give  rise  to  the  same  direction  of  flow  at  every  point  on  a  hemisphere.  The  only  possible 
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configurations  of  motions  that  may  be  contradictory,  provided  the  surfaces  in  view  satisfy 
the  constraints  described  in  Section  3,  are: 

a:  three  or  more  motions  such  that  the  corresponding  zero  motion  contours  gu  =  0  inter¬ 
sect  in  the  same  point  p^,  with  utr.  (Pij)  x  ^trj(Pij)  7^  0  for  z  j. 

b:  three  or  more  motions,  such  that  all  corresponding  zero  motion  contours  ga  =  0  are 
tangential  at  the  same  point  ri2,  which  as  described  in  Appendix  E,  can  occur  only  if 
—  Q  for  some  constant  C. 

tanZ.(Ct’i,ri2) 


Appendix  B  Zero  motion  contours 


Let  us  consider  the  following  question:  What  is  the  locus  of  points  where  the  flow  due  to 
the  given  rigid  motion  can  possibly  be  zero?  As  in  [4]  we  can  show  that  such  points  are 
constrained  to  lie  on  a  second  order  curve  on  the  sphere. 

The  flow  at  point  r  can  be  zero  only  if  the  rotational  and  translational  components  at 
r  are  parallel  to  each  other.  Let  t  and  w  be  translational  and  rotational  velocity  of  the 
observer.  Then  the  flow  at  point  r  can  be  zero  only  if 

(r  X  (t  X  r))  X  (a?  X  r)  =  0  (34) 


By  simple  vector  manipulation,  from  (34)  we  obtain 


((w  X  r)  •  (t  X  r))r  =  0 

Since  r  7^  0,  the  flow  at  point  r  can  be  zero  only  if 

(a;  X  r)  •  (t  X  r)  =  0 

or  07  •  t  —  (w  •  r)(t  •  r)  =  0 


(35) 


(36) 


Equation  (36)  describes  a  second  order  curve  on  the  sphere,  which  we  will  call  the  zero 
motion  contour  of  the  rigid  motion  (t,o7).  The  zero- motion  contour  consists  of  two  closed 
curves  on  the  sphere.  As  shown  in  Figure  14,  if  (07  -t)  >  0,  one  of  the  curves  contains  to  =  -jli 
and  070  =  and  one  contains  —to  and  — cjq;  if  (u7  •  t)  =  0  the  two  curves  become  great 
circles,  one  orthogonal  to  t,  the  other  orthogonal  to  07;  if  (u7  •  t)  <  0  one  of  the  two  curves 
passes  through  to  and  — a7o  and  the  other  through  —to  and  079. 
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Figure  14:  The  zero  motion  contour  (the  locus  of  points  r  where  r  could  be  zero)  consists 
of  two  closed  curves  on  the  sphere.  Three  possible  configurations  are  (a)  (w  •  t)  >  0,  (b) 
(w  •  t)  =  0,  and  (c)  (a?  •  t)  <  0. 

Appendix  C  Zero  motion  contours  are  not  tangent 

To  show  that  zero  motion  contours  gw  and  gw  are  not  tangent  at  ri  (see  Section  4.1,  part  A), 
let  us  compute  tangent  lines  to  the  contours  at  ri.  Let  the  direction  of  the  line  tangent  to 
^11  at  ri  be  Ui.  The  line  lies  in  the  plane  tangent  to  the  sphere,  so 

Ui  •  ri  =  0  (37) 

Directional  derivative  of  511  along  Ui  must  be  zero.  Let  r^  =  ri  +  Ui£.  Then 

=  (oji  X  ui)  •  (ti  X  ri)  +  (wi  X  ri)  •  (ti  x  ui)  +  x  ui)  •  (ti  x  Ui)  (38) 

as 

Since  (37)  holds,  Ui  must  also  satisfy 


-  (uJi  •  ri)(ux  •  ti)  -  (wi  •  Ui)(ri  •  ti)  =  0 

(39) 

Thus  we  obtain 

ui  •  ((wi  •  ri)ti  +  (ti  •  ri)u;i)  =  0 

(40) 

We  can  compute  the  tangent  direction  from  (40)  and  (37)  as 

ui  =  ri  X  ((wi  •  ri)  ti  +  (ti  •  ri)  wi) 

(41) 
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Similarly,  the  direction  tangent  to  gi2  is 


U2  =  ri  X  ((u>i  •  ri)t2  +  (t2  •  ri)a;i) 

Since  point  ri  lies  on  geodesic  n,  we  get 

U1XU2  =  0  + (wi  •ri)(t2  •ri)((ri  X  ti)  X  (ri  X  wi))+ 

(ti  •  ri)(a;i  •  ri)((ri  x  wi)  x  (ri  x  t2))  +  0  = 

=  (wi  •  ri)((t2  •  ri)[t<^i  ri  ti]  rj  +  (ti  •  ri)[t2  ri  oji]  ri) 

=  (a>i  •  ri)((ti  X  t2)  •  (ri  x  (wi  x  ri)))  ri 

Also 

a>i  •  ri  =  llwi  X  (ti  X  t2)|p  >  0 

and 

(ti  X  t2)  ■  (ri  X  (wi  X  ri))  =  [tit2UJi]  ||a;i  x  (ti  x  t2)||^  ||ti  x  t2f  ^  0 
So  Ui  X  U2  is  not  zero,  and  the  two  zero  motion  contours  cross  at  point  ri. 

Appendix  D  Zero  motion  contour  crossing  the  border  of  the  image 

Let  the  half  sphere  image  be  bounded  by  equator  let  no  be  a  unit  vector  normal  to  the 
plane  of  q.  We  would  like  to  know  whether  the  zero  motion  contour  of  motion  (t,  uj)  intersects 
equator  q. 

Let  us  choose  a  Cartesian  coordinate  system  such  that  no  =  [0,0,1].  Let  t  = 
and  w  =  [uJx,ojy,LOz]-  Points  on  equator  q  can  be  written  as  [cos^,  sin^,  0].  Thus  the  zero 
motion  contour  (u;  x  r)  •  (t  x  r)  =  0  intersects  q  if  equation 

{iOxix  +  '^ztz)  sin^  (f)  ~  {^yix  +  l-^x'ty)  siu  (f)  cos  6  +  {ijJyty  +  OJztz)  cos^  ^  =  0  (46) 

has  a  solution. 

Writing  A  =  tan  we  obtain  a  quadratic  equation 

T  ^ ij^yix  T  ^x^y')^  d”  {y^y^y  T  ^z'^z'}  —  0  (4T) 


(42) 

(43) 

(44) 

(45) 
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This  equation  has  a  real  solution  if 


I  -  {^Oytx  "h  ^X^y)  H"  U^gtz){y^yty  4"  ^z^z)  ^  0 

After  some  manipulation,  we  obtain 

I  —  {iOytx  i^xty')'^  ^{^z^z'){^xtx  4“  ^y^y  4"  ^z^z) 

—  [twno]^  -  4(a>  •  no)(t  •  no)(t<;  •  t)  >  0 

Appendix  E  Intersections  of  zero  motion  contours 

Let  (ti,Wi)  and  (t2,  W2)  be  two  ambiguous  motions.  Let  us  investigate  possible  intersection 
points  of  the  zero  motion  contours  of  the  two  motions. 

Since  ambiguity  is  possible  only  when  Wi  xa;2  7^  0,  ti  xt2  7^  0,  and  (uJi  xti;2)‘(ti  xt2)  =  0, 
we  can  choose  a  Cartesian  coordinate  system  such  that 

X  =  («i  X  a;2)/|wi  X  a;2| 

y  =  (ti  X  t2)/|ti  X  t2|  (50) 

Z  =  XxY 

In  this  coordinate  system,  we  can  write  ti  =  [£/i,0,  Wi],  t2  =  [I72?0,  W2],  u)\  =  [0,,5i,7i], 
and  UJ2  =  [0,/?2,72]- 

Clearly,  both  zero  motion  contours  pass  through  point  [0,0, 1].  We  would  like  to  know 
whether  this  is  the  only  intersection  point. 

If  Wi  were  zero,  we  would  have  ti  •  =  0;  thus  the  zero  motion  contour  would  be 

degenerate.  This  is  not  possible  if  the  motions  are  ambiguous.  Thus  Wi  7^  0,  and  similarly 
14^2,  7i,  and  72  non-zero. 

Since  the  zero  motion  contour  does  not  depend  on  the  size  (and  direction)  of  vectors  t 
and  oj,  we  can  re-scale  vectors  ti,  Wi,  multiplying  by  Aj  7^  0  such  that 

ti  =  Ai  ti  =  [[7^1,0, 1] 

t2  =  A2t2  =  [6^2,0,!] 

Wi  =  AsUJi  =  [0, /3i,l] 

0^2  =  A4W2  =  [05^25!] 


(48) 


(49) 
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Let  us  consider  point  r  =  [x^y^z]  such  that  (w  x  r)  •  (t  x  r)  =  0.  li  z  ^  0,  point 
r/z  =  [a;/z,y/z,  1]  also  satisfies  the  equation.  Thus  it  is  enough  to  consider  two  possible 
sets  of  points;  points  of  the  form  r  =  [x,y,  1]  (i.e.,  points  lying  in  the  plane  tangent  to  the 
sphere  at  [0,0,1]),  and  points  r  =  [a:,y,0]  (these  points  correspond  to  points  at  infinity  on 
the  tangent  plane). 

A:  To  obtain  the  possible  intersection  points  r  =  [a;,?/,  1]  we  express  the  zero  motion 
contours  as 

(wi  X  r)  •  (ti  X  r)  =  -  xU-^,  -  y^i  -  xyUi^i  -  0  (52) 

(^2  X  r)  •  (t2  xr)  =  x^  +  y‘^  -  xU2  -  y/§2  -  xyU2^2  =  0  (53) 

We  can  compute  y  from  the  difference  of  the  two  equations  as 

xjU,  -  U,) 

^  x{ij2h-Uih)^h-h 

Substituting  (54)  into  (52),  we  obtain  a  polynomial  equation  of  degree  4  in  x.  One 
solution  is  a:  =  0  (both  zero  motion  contours  pass  through  point  [0,0,1]).  The  remaining 
equation  of  degree  3  has  at  least  one  real  solution. 

If  01  —  ^2){U02  —  U2^i)  =  0,  we  obtain  another  solution  x  =  0.  Otherwise  the  two 
contours  intersect  in  two  different  points.  Since  uji  x  u;2  ^  0,  we  know  —  ^2  ^  0.  Thus 

the  two  zero  motion  contours  are  tangent  only  if  Ui(32  =  U2$i-  If  this  is  the  case,  we  obtain 

an  equation  of  degree  2  in  a.  Its  discriminant  is 

d  =  PlU^{U2  +  UifiWif  -  4)  (55) 

Thus  if  \U2/3i\  =  \Ui^2\  >  2,  the  two  zero  motion  contours  are  tangent  at  [0,0,1],  but 
intersect  at  two  other  points. 

B:  Now  we  compute  intersection  points  r  =  [a;,z/,0].  We  can  assume  =  1.  Then 

we  obtain  equations  xyUi  =  1  and  xy  U2  ^2  =  T  so  such  intersection  point  exists  only  if 

U0X  =  U2k- 

In  the  previous  part  we  have  shown  that  if  U2$i  7^  Ui$2^  the  two  zero  motion  contours 
have  more  than  one  intersection  point.  So  it  is  enough  to  check  the  tangential  case  here. 
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If  the  two  contours  are  tangent  at  [0,0, 1],  from  =  U2P2  —  Ui^2  we  obtain 

|t^i|  =  \U2V  Since  ti  x  t2  ^  0,  this  is  possible  only  if  Ui  =  -IJ2  and  ^1  =  -^2- 
Writing  x  =  cos  (j)  and  y  =  sin  we  obtain  equation 

2  2 


sin(2^)  = 


(56) 


so  again  there  is  an  intersection  point  ii  \Ui^2\  > 

Therefore  the  two  zero  motion  contours  have  only  one  intersection  point  if  =  C/2/^2, 

and  \Ui^2\  =  \U2^i\  <  2.  If  we  denote  the  intersection  point  of  5^11  and  g22  as  ri2,  this  can 
be  written  as 

(57) 


tanZ(ti,ri2)  tanZ(t2,ri2) 


tanZ(wi,ri2)  tan  Z(a;2,  ri2) 
where  Z(-,-)  denotes  the  angle  between  two  vectors. 

This  relationship  can  also  be  expressed  as 


[a7iU?2tl]  [tit2th2]  =  [a;ia>2t2][tlt2d?i] 

Furthermore  from  (55)  we  obtain  the  constraint 

[I’lf2^2]|  11^1  ^  ^2!!  Il^l  ^  I’2|| 

2  |(((a;i  X  UJ2)  X  (ti  x  t2))  •  ti)  (((wi  x  ^2)  x  (ti  x  tg))  •  u;2)| 

Appendix  F  Antipodal  contradictory  points 

Here  for  the  purpose  of  providing  a  description  of  the  areas  where  two  motion  fields  are 
ambiguous,  the  conditions  are  developed  for  point  r  and  its  antipodal  point  -r  both  to  be 
contradictory. 

Clearly,  if  one  of  the  CP-conditions  holds  for  r,  it  cannot  be  true  for  — r.  So  if  both  r 
and  —  r  are  contradictory,  two  of  the  conditions  must  hold  at  r  and  the  other  two  at  —  r. 

If  (18)  and  (19)  hold  at  r  and  (20),  (21)  at  -r,  we  get 

sgn(/i(r))  =  sgn(5i2)  =  -sgn(^22)  =  sgn(^2i)  =  -sgn(/u,(r)) 

at  r  and 

-sgn(/f(r))  =  -sgn(5-2i)  =  sgn(^n)  =  -sgn(^i2)  =  sgn(/>^(r)) 


(58) 

(59) 
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at  — r.  Thus 


Sgn(^ii)  =  sgii(^22)  =  -Sgll(^12)  =  -sgn(52i) 
sgn(/f(r))  ^  sgn(/^(r)) 

If  (18)  and  (20)  hold  at  r  and  (19),  (21)  at  -r,  we  get  sgn(5(i2)  =  -sgn(5'22)  at  r  and 
sgn(5i2)  =  sgn(5r22)  at  — r,  so  this  case  cannot  occur. 

If  (18)  and  (21)  hold  at  r  and  (19),  (20)  at  — r,  we  get 

sgn(/t(r))  =  sgn(<?i2)  =  -sgn(^22)  =  sgn(/^(r))  =  -sgn(^ii) 

at  r  and 

-sgn(/^(r))  =  -sgn(^2i)  =  sgn(522)  =  -sgn(/t(r))  =  sgn(^ii) 
at  — r.  So  we  get 

sgn(5'n)  =  sgn(5r22)  =  -sgn(5'i2)  =  -sgn(p2i) 
sgn(/t(r))  =  sgn(/,,(r)) 

Thus  point  r  and  point  — r  are  both  contradictory  if  and  only  if 

sgn(^ii(r))  =  sgn(5r22(r))  =  -sgn(5i2(r))  =  -sgn(^2i(r))  (60) 
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